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Abstract 

The technique of applying form- invariant, spatial coordinate transformations of 
Maxwell's equations can facilitate the design of structures with unique electromag- 
netic or optical functionality. Here, we illustrate the transformation-optical approach 
in the designs of a square electromagnetic cloak and an omni-directional electro- 
magnetic field concentrator. The transformation equations are described and the 
functionality of the devices is numerically confirmed by two-dimensional finite ele- 
ment simulations. The two devices presented demonstrate that the transformation 
optic approach leads to the specification of complex, anisotropic and inhomogeneous 
materials with well directed and distinct electromagnetic behavior. 
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1 Introduction 



In a theoretical study, Pendry et al. reported a general method for the de- 
sign of electromagnetic materials based on form-invariant transformations of 



Preprint submitted to Elsevier 



1 February 2008 



Maxwell's equations [1]. In that paper, the methodology of transformation 
optics was applied to find the specification for an electromagnetic cloak-a 
complex material capable of rendering objects within its interior invisible to 
detection. Although just one example of the many intriguing structures pos- 
sible using the transformation optical approach, the proposed cloak design 
generated enormous interest in its own right. An approximation to the invis- 
ibility cloak based on metamaterials was subsequently realized by Schurig et 
al., who demonstrated the cloaking mechanism in microwave experiments [2]. 
The transformation optical approach to invisibility is quite general, differing 
in scope from prior related work. Indeed, methods of reducing the electromag- 
netic scattering of objects at radar frequencies have long been a subject of 
intense research [3,4,5]. On the nanoscale, techniques have also been suggested 
to reduce the scattering of one or more multipole components of size-limited 
objects using tailored negative index or negative permittivity coatings [6,7]. 
More recently, a mathematically rigorous proof of an invisibility structure 
based on active devices was reported [8]. 

Transformation optics provides for a conceptually simple approach to the de- 
sign of complex electromagnetic structures: one imagines warping space to 
achieve the desired electromagnetic functionality. The trajectories of electro- 
magnetic waves passing through a region of warped space must conform to 
the local metric, and this provides an alternative (though conceptual) means 
to control and manipulate electromagnetic fields. Once the desired design 
is determined, the coordinate transformation and its Jacobi matrix deter- 
mine the transformation of Maxwell's equations and the constitutive rela- 
tions. The result provides the specification for an electromagnetic structure 
that is complex-being inhomogeneous and anisotropic-but realizable for exam- 
ple through artificially structured metamaterials. Indeed, because the fields in 
a volume bounding a transformation optical structure are identical to those 
that would exist where the structure is replaced by free space, anisotropy is 
necessary to circumvent uniqueness constraints [8]. 

If the coordinate transformation can be realized exactly in the constitutive 
parameters, all aspects of wave propagation will be transformed by the struc- 
ture, including the near-fields. Adding constraints to the materials reduces 
the ultimate performance of the structure, but nevertheless can still allow for 
interesting and novel structures. Leonhardt, for example, has shown that if 
the materials are restricted to be isotropic, an approximate cloak can be con- 
structed that is valid in the geometrical optic limit [9]. Likewise, constraints 
were employed for the metamaterial cloak utilized by Schurig et al. to ease 
the metamaterial design and fabrication, resulting in a structure that pro- 
duced significant reflection yet still demonstrated the cloaking mechanism for 
transmitted waves [2]. 

Since the concept of transformation optics was introduced, there have been a 
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growing number of subsequent reports applying the method to a variety of elec- 
tromagnetic, acoustic and elasto-mechanical structures [10,11,12,13,14,15,16,17] 
Full wave simulations have helped to confirm the expected behavior and have 
provided a platform to explore systematically the effects of absorption, imper- 
fections and other constraints that are inherent to fabricated realizations of 
the transformation optical structures [18,19]. 

In this paper we present two examples that demonstrate the general applica- 
bility of form-invariant coordinate transformations for the design of complex, 
inhomogeneous and anisotropic electromagnetic materials with well-defined 
functionality. For the first example, we derive the electromagnetic constitutive 
parameters corresponding to a two-dimensional electromagnetic cloak having 
square cross-section. The square shape has been chosen to illustrate the na- 
ture of the transformation-and the resulting design-for a structure that lacks 
rotational symmetry in the plane. In contrast to the cylindrical cloaks with 
circular cross-section previously presented, the square cloak design results in a 
non-orthogonal transformation producing a more complicated specification for 
the spatially dependent permittivity and permeability tensors. The method to 
design this structure, however, can be applied to the design of structures with 
arbitrary shape. 

For the second example, we derive the material properties of an electromag- 
netic field concentrator by the same approach. The purpose of the cylindrical 
concentrator is to focus incident electromagnetic waves with wave vectors per- 
pendicular to the cylinder axis, enhancing the electromagnetic energy density 
of incident waves in a given area. This example illustrates the strength of the 
transformation-optical approach for designing devices other than cloaks. 



2 Transformation Equations 

In this section, the formulas describing the spatial coordinate transforma- 
tions and the calculation of the resulting material parameters, i. e. the electric 
permittivity tensor and the magnetic permeability tensor, are derived. The 
methodology used to compute the electromagnetic material properties is sim- 
ilar to the one reported in [20] . 

For convenience, we denote Maxwell's equations in the Minkowski form [21] 



where the square brackets express an alternation among the indices [22] and 
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the skew-symmetric covariant and contravariant tensors F\ v and G uX and the 
contravariant vector j x possess the identifications 
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where E is the electric field, B is the magnetic induction, D is the electric 
displacement, H is the magnetic field, p is the volume charge density, j is 
the current density and the indices denote their spatial components. In this 
notation, the coordinate vector in four-space is x a = (xo = t,xi,x 2 ,x 3 ) T with 
the vacuum speed of light c set to unity. 



For a linear medium, the constitutive relation can be written as 

1 
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where the tensor ■y XvuK contains the complete information about the electro- 
magnetic material properties. 

The Minkowski equations (1) and (2) and the constitutive relation (6) are 
form-invariant for arbitrary continuous space-time transformations of the form 



x a \x a ) = A^x a 



(7) 



where A: 



dx° 



dx a 



are the elements of the Jacobian transformation matrix and 



the primed indices denote the space-time coordinates of the vector x in the 
transformed space. 
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Considering the transformation (7), the Minkowski equations (1) and (2) 
transform as 



d[n'F\' v >] = A K K ,A x x ,A»,d [K F Xv] = (8) 
d u ,G u ' x ' = [det(A x ')]- 1 A*d 1/ G l/X 

= [det(A^A x :f (9) 

= / (10) 
and the constitutive relation transforms as 
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with 
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where det(X) indicates the determinant of a tensor X. 

It is obvious, that the form-invariance of the Minkowski equations (1) and 
(2) and the constitutive relation (6) also holds for transformations which only 
address the space coordinates, as the space manifold is a sub-manifold of 
the space-time manifold. In the following we will restrict ourselves to time- 
independent, spatial coordinate transformations. Under this restriction, the 
constitutive parameters, i. e. the tensors of the permittivity and the perme- 
ability of a linear, anisotropic, non- dispersive, non-bianisotropic and locally 
interacting medium can be written in a more accessible form as 



j'l' = [det(JlZ)]- 1 A*A>' ( # (14) 

/Mdet^r^y (15) 

The relations (7), (14) and (15) form the underlying equations for the cal- 
culation of the electromagnetic material parameters used in the design of a 
square-shaped cloak and a concentrator for electromagnetic fields. 

For all the transformations considered in the next sections, the mathematical 
starting point is three-dimensional, euclidian space expressed in a cartesian 
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coordinate system x l = (x\, x 2 , x 3 ) T . From the physical point of view, the space 
is considered to be medium-free (vacuum) and isotropic. Thus, the permittivity 
tensor e u and the permeability tensor of the original space can be expressed 
in the form 



e lj =e 5 i] (16) 
// J =/io^' (17) 

with (5 lj = 1) for (i = j) and (5 lj = 0) elsewhere. 



2.1 Square Cloak 



The coordinate transformation equations for the electromagnetic design of a 
square-shaped cloak with a sidelength 2s\ of the inner square and a sidelength 
2s 2 of the outer square (see fig. la) are expressed by 



x' 1 (x 1 ,x 2 ,x 3 )=xi— — — + s 1 (18) 

x' 2 (x 1 ,x 2 ,x 3 )=x 2 (— — — + — ) (19) 
V s 2 xiJ 

x' 3 (x 1 ,x 2 ,x 3 )=x 3 (20) 
with the Jacobi matrix and its determinant 
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for (0 < x\ < s 2 ), (—52 < x 2 < s 2 ), \x 2 \ < \xi\ and (|xs| < oo). It should be 
noted, that, by the foregoing definitions, the transformation equations (18)- 
(20) are only defined for the green shadowed area in fig. la and that the 
transformation is continuous at the boundary of the transformed domain. The 
corresponding transformation formulas for the upper, left, and lower domain of 
the square cloak can be readily obtained by applying rotation operators with 
rotation angles of ir/2, ir and 3n/2 around the z-axis to equations (18)-(20). 

As can be seen from fig. lb, the transformation expands the space within 
the inner square at the expense of a compression of space between the inner 
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Fig. 1. Spatial coordinate transformation for the design of a square-shaped cloak 
(a) original space, s\: half sidelength of the inner square, S2- half sidelength of the 
outer square. The transformation equations (18)-(20) are only valid in the green 
shadowed region (b) transformed space 



and outer square. Inserting equations (21- 22), (16-17) into (14) and (15), 
applying the reverse transformations x t (x l ) of (18)-(20) and dividing (14) by 
Eq and (15) by [1q yields the relative permittivity and the relative permeability 
tensors (e r ) 1 J and (/i r ) lJ , expressed in the coordinates x l of the transformed 
space, as 
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Fig. 2. (a) Rotation angle of the eigenbasis vectors (x, y) in dependence on their 
spatial location (x\,X2) within the cloaking material calculated for the green shad- 
owed region in fig. la. (b-d) Spatial dependence of the material parameters of the 
square cloak medium for a TE-wave polarized in z-direction, expressed in their lo- 
cal eigenbasis (x,y,z), (b) permeability fj, x in x-direction, (c) permeability fi y in 
y-direction, (d) permittivity in e z in z-direction 
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with 



.S'2 , Xo ( Sj I . . 

a:= 7^7^ h '- = W aSuC:=a \ (24) 

Due to the natural invariance of the Minkowski equations (as discussed in 
[20]), the permittivity and permeability tensors can also be interpreted as the 
material properties of a medium described in the coordinate system of the 
original space by substituting the primed indices by unprimed indices ("ma- 
terial interpretation"). Again it should be noted, that the material properties 
(23) are only valid in the green shadowed region of fig. la and that, due to 
the symmetry of the cloak, the material properties of the other cloak domains 
can be readily obtained by rotating the tensors in (23) by ir/2, n and 37r/2, 
respectively. Furthermore, the relative permittivity and permeability tensors 
are non-diagonal, which is a direct consequence of the non-conformality of the 
transformation (18)-(20). However, in terms of fabricating such a material, it 
is desirable to have the material parameters denoted in their eigenbasis, where 
the permittivity and permeability tensors are diagonal. Due to the symmetry 
of the tensors e v and fi^ , an eigenbasis solution always exists and one obtains 
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with 



A: =a 2 + b 2 + c 2 

B := a 2 [a 2 + 2{b 2 - c 2 )) + b 2 {b 2 + 2c 2 ) + c 4 



(26) 
(27) 



Notice, that the primes were omitted to express the electromagnetic param- 
eters in the material interpretation. Due to the spatial dependence of the 
elements of the permittivity and permeability tensors, the orientation of the 
basis vectors (x,y) of the eigenbasis depends on the spatial location (xi,x 2 ) 
within the cloak material. This is illustrated in fig. 2a for a square cloak with 
Si = 3 cm and s 2 = 6 cm. The graph shows the rotation angle of the eigen- 
basis vectors (x,y) in dependence on the location x\ at different positions x 2 . 
Again, the physical quantities are calculated for the green shadowed region 
in fig. la. The rotation angle of the eigenbasis with refer to the coordinate 
system (xi,x 2 ) varies within a range from 1.3 to 1.8 rad. In order to fabricate 
such a medium as a metamaterial, the principle axes of the unit cells have to 
be individually aligned along the basis vectors (x, y, z) of the eigenbasis. 

Figures 2b-d show the values of the relative permeabilities fi x := fij 1 and 
H y := /i 22 and the relative permittivity e z := e 33 in dependence on the location 
within the cloak. These three physical quantities deliver a full description of 
the propagation behavior of an electromagnetic wave with a linear polarization 
vector of the electric field oriented along the z-direction. The depicted area 
corresponds to the green shadowed region in fig. la. As opposed to fi y and 
e 2 , the value of \i x diverges at the boundary of the inner square of the cloak. 
However, at a distance of about 1.6 mm from the inner boundary of the square 
cloak, fi x already approaches finite values below 35 as determined along a 
straight intersection line parallel to the y-axis at x — 3.16 cm. Please note, 
that fig. 2b only displays the permeability \i x for x > 3.16 cm and thus does not 
show the divergence of the relative permeability towards the inner boundary of 
the cloak. However, assuming a typical unit cell size of 3.3 mm for a fabricated 
effective medium at a working frequency of 8.5 GHz, the value of \i x at the 
midpoint of the unit cells at the inner boundary of the cloak is between 19 and 
35, so that an implementation of such a material is still possible. Although the 
effective fi x is inaccurate in the vicinity of the boundary of the inner square, it 
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can be shown, that the performance of the implemented device is not affected 
by this fact, which is out of the scope of this paper. 



2.2 Cylindrical Concentrator 



Due to its cylindrical symmetry, it is convenient to describe the transformation 
equations in a cylindrical coordinate system. In this context, it is necessary to 
consider the transformation from cartesian to cylindrical coordinates for an 
isotropic medium with permittivity e and permeability ix. With the transfor- 
mation 



r'(x 1 ,x 2 ,x 3 ) -- 
(f)'(xi,X2, x 3 ) = arctan ^ 
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and equations (14) and (15) one obtains 
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At this point, the reader should be aware, that the metric tensor of the trans- 
formed space is 
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so that (31) necessarily still represents the material properties of an isotropic 
medium. 

The transformation equations for the optical design of the cylindrical concen- 
trator are denoted as 
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Fig. 3. Visualization of the space transformation expressed in (33)-(35). Space is 
compressed into a cylindrical region with radius R\ = 2 cm (black inner circle) 
at the expense of an expansion of space between R\ and R 3 = 6 cm (black outer 
circle). The displayed intermediate circle is located at R2 = 4 cm. 
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with the corresponding Jacobi tensors and determinants 
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< r' < R 2 
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The space transformation is visualized in fig. 3. Space is compressed into 
a cylindrical region with radius Rl at the expense of an expansion of space 
between Rl and R3. The transformation is continuous to free space at R3. 
Inserting (31), (36) and (37) into (14-15) and renormalizing (31) by requir- 
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Fig. 4. Radial dependence of the material properties of a cylindric concentrator with 
an inner radius R\ = 2 cm and an outer radius R2 = 6 cm. (a) radial permeability 
component ji r , (b) azimuthal permeability component /j,^, (c) z-component e z of the 
permittivity 

describe the relative material properties of free space as rf^' = 8 1 ^' in cylin- 
drical coordinates, one obtains with help of the inverse transformations of 
(33-35) the relative permittivity and permeability tensors, expressed in the 
coordinates (r", <fi", z") as 
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with 
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(3* 



(39) 
(40) 



e := R 2 — Ri, f := R3 — R2, h := R 3 — Ri 

Again, in the material interpretation, (38) represents the material properties 
of the cylindric concentrator in the original space (r f , <f)', z'). 



,2'2' 



Fig. 4 shows the radial and azimuthal components fi r := fil and \i$:= \i 
and the z-component e z := 3 in dependence on the radial position within 
the concentrator, assuming an interaction with an electromagnetic wave with 
a polarization of the electric field parallel to the z-direction. The experimental 
implementation of such a material requires independent control of the local 
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Fig. 5. Computational domain terminated by perfectly matched layers with func- 
tionality in the x-, y- and x-y-direction. The TE plane wave is excited by a current 
sheet. 

values of all three parameters. The metamaterial design and fabrication is part 
of current research. 



3 Simulation Results and Discussion 

For the full wave electromagnetic simulations a two-dimensional finite-element 
solver of the Comsol Multiphysics software package was used. The computa- 
tional domain and its boundaries are shown in fig. 5. A transverse-electric (TE) 
plane- wave was excited by a current sheet. The computational domain was ter- 
minated by perfectly matched layers (PMLs). The current density distribution 
was chosen to exponentially decrease towards the borders of the sheet in the 
y-direction to prevent interactions with the PMLs. The calculations did not 
include absorption in the object. The finite-element solver required all phys- 
ical quantities to be described in cartesian coordinates. For all simulations, 
a stationary solver was used. The solver allowed to conveniently implement 
the functional dependence of the permittivity and permeability tensors of the 
simulated material into the model and thus to accurately describe and predict 
the electromagnetic behavior of the designed medium. 

3.1 Square Cloak 

Fig. 6 shows the results of the two-dimensional full-wave simulations of a 
square-shaped cloak. The color map depicts the spatial distribution of the 
real part of the transverse-electric phasor oriented along the z-direction. In 
addition, the direction of the power flow is indicated by the grey lines. The 
frequency of the TE wave is 8.5 GHz. The sidelengths of the inner and outer 
square of the cloak are 6 cm and 12 cm, respectively. 
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Fig. 6. Electric field distribution in the interior and exterior region of the square 
cloak. The direction of the power flow is in the positive x-direction. The wave is 
smoothly bent around the corners of the square cloak, (a) Phase fronts parallel to 
one side of the cloak, (b) Cloak rotated by tt/8 with respect to the phase fronts of 
the incoming wave 



In fig. 6a, the phase fronts of the impinging wave are parallely aligned to one 
of the sides of the cloak. As can be seen, the wave is smoothly bent around 
the cloaked area and the phase fronts are completely restored when the wave 
exits the cloak material. The inhomogeneity and the anisotropy of the cloak 
medium are evident as the direction of the power flow and the phase front 
normal are not parallel and the angle between the directions changes locally. 
In fig. 6b, the cloak is rotated by an angle of ir/8 with respect to the phase 
fronts of the incoming wave. In this configuration, the phase fronts are no 
longer parallel to any side of the square cloak. As before, the phase fronts 
are completely restored after propagation through the cloak material and the 
inner square is not sensed by the wave. In both cases, the wave impedances 
of the cloak medium and free space are exactly matched and the device is 
therefore reflectionless. 



The square-shaped cloak is an example of a cloak with reduced symmetry 
in comparison to a cylindric cloak. In addition, the square cloak possesses 
sharp corners. The simulations clearly show, that the transformation-optical 
cloak design is not restricted by cylindric symmetry requirements. In principle, 
cloaks of arbitrary shape can be designed by use of the transformation-optical 
approach. 
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Fig. 7. (a) Electric field distribution in the interior and exterior region of the cylin- 
dric concentrator. The direction of the power flow is in the positive x-direction, 
indicated by the grey lines, (b) Normalized power flow distribution. The power flow 
is enhanced within the region with radius Rl=2 cm by a factor of 2. Much stronger 
power flow enhancements can be achieved by increasing the ratio R2/R1. 

3.2 Concentrator 



In section 2.2, the material properties of a concentrator were described in a 
cylindrical coordinate system. In cartesian coordinates, the relative permittiv- 
ity and permeability tensors can be obtained from (38) by use of the general 
transformation 



(e'-e 2 )^ ^(e'xj+e'xi) o (41) 
e% 

with £ = e r , fi r . The variables e r and fi r as functions of the space variables 
(x := x%,y := X2, z := 23) are then directly assigned to the concentrator 
domains. 
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Fig. 7a displays the real part of the phasor of the electric field for a z-polarized 
TE wave. The grey lines represent the direction of the power flow. The free- 
space frequency of the TE wave is 8.5 GHz. The outer radius of the concentra- 
tor is #3 = 6 cm. As can be seen, the fraction of the plane-wave extending in 
the y-direction from — R 2 = —4 cm to R2 = 4 cm is completely focussed by the 
concentrator into the region with radius R\ = 2 cm. Additionally, the fields 
within the intervals [— R 3 , — R2) and (R^jRs] in the y-direction are focussed 
to an area with a radius lying in the interval (Ri, R 3 }. 
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Fig. 7b illustrates the normalized intensity distribution of the TE wave. It is 
obvious, that the field intensities are strongly enhanced in the inner region of 
radius R\ within the concentrator material. The intensity enhancement factor 
for the chosen structure, computed as the ratio between the maximal values 
of the field intensities outside the circular region with radius R3 and inside 
the concentrator region with radius R\, is 2. Significantly stronger enhance- 
ments can be achieved by increasing the ratio i?2/-Ri- As can be seen, the 
enhancement theoretically diverges to infinity as Ri goes to zero. 

Due to the rotational symmetry around the axis perpendicular to the x-y- 
plane, the concentrator focusses waves impinging from arbitrary directions. 
The concentrator is reflectionless due to inherent impedance matching in 
the transformation-optical design method. Although metamaterials, which are 
necessary to implement the material properties of a concentrator, inherently 
suffer from losses, we think, that the concentrator can play an important role 
in the harnessing of light in solar cells or similar devices, where high field 
intensities are required. 



4 Conclusion 

In conclusion, we have presented the material design of a square-shaped cloak 
and an electromagnetic field concentrator based on form-invariant transfor- 
mations of Maxwell's equations. The electromagnetic behavior of the devices 
was simulated by use of a two-dimensional finite element solver. In contrary 
to previous publications, the simulated cloaking device did not possess a 
cylinder-symmetry. The proposed electromagnetic field concentrator proved to 
be well suited for the confinement of electromagnetic energy of waves imping- 
ing from arbitrary directions. The two demonstrated optical devices exemplify 
the strength of the general methodology of form-invariant coordinate trans- 
formations of Maxwell's equations for the design of electromagnetic materials 
with a well-defined functionality. The technique allows to chose from an infi- 
nite set of allowed transformations and thus provides a powerful tool for the 
conception of optical elements with previously unachievable electromagnetic 
behavior. 

D. Schurig wishes to acknowledge support from the IC Postdoctoral Research 
Fellowship program. 
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